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Abstract 

The coalescence of three levels has particular attractive features. Even though it may be difficult 
to realise such event in the laboratory (three additional real parameters must be adjusted), to take 
up the challenge seems worthwhile. In the same way as the chiral behaviour of a usual EP can give 
a direction on a line, the state vectors in the vicinity of an EP3 provide an orientation in the plane. 
The distinction between left and right handedness depends on the distribution of the widths of the 
three levels in the vicinity of the point of coalescence. 



1 



INTRODUCTION 



There is substantial literature during the past decade relating to exceptional points [lj], 
i.e. points where two eigenvalues of an operator coalesce giving rise to a square root singular- 
ity in the spectrum Q . These singularities are generic and are thus encountered in virtually 
all physical problems associated with eigenvalues. They have been discussed in mechanical 



problems [3], in optics Ml, for bound states [5| and resonances |6|, |7[ in quantum mechanics 
and in atomic physics [8|. The mutual influence of neighbouring exceptional points upon 
the phase behaviour of the associated wave functions is dealt with in 9|. Exceptional points 
also play a crucial role in quantum phase transitions [h]]. Chiral behaviour of the eigen- 
function as well as effects of time reversal symmetry breaking is discussed in [3, 
Experimental manifestations including chiral behaviour have been achieved with microwave 
cavities ||14l and coupled oscillators in electronic circuits 15]. A recent more mathemat- 

n 

ical expose investigates exceptional points in the context of projective Hilbert spaces [16( 
while a connection to PT symmetric Hamilton operators 17| is found in [18]. In fact, it 
is established in [l^] that, for a pseudohermitian VT symmetric Hamiltonian, the onset of 
spontaneous symmetry breaking by the wave function happens just at an exceptional point. 

In many of the papers quoted above the notation EP was used denoting the generic 
exceptional point where two levels coalesce at a square root branch point. To distinguish 
it from the non-generic coalescence of three levels being the subject of the present paper 
we'll denote the generic EP by EP2 and we denote by EP3 the specific situation where three 
levels coalesce. We note that EPs of higher order have been implicitly encountered as the 



0,6 



coalescence of two or more EP2 in a recent investigation of a complex WKB analysis 



20). 



THREE LEVELS COALESCING 



The situation of three or more levels coalescing does not seem to have been investigated 
in great detail, the reason being that there are too many parameters needed to enforce 
such higher order coalescence. In fact, while two real parameters (one complex parameter) 
suffice to invoke the coalescence of two levels, for N levels coalescing (N 2 + N — 2)/2 real 
parameters are needed considering complex symmetric matrices. For N = 3 it means that 
three additional real parameters have to be chosen judiciously to invoke the coalescence 
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of three levels in the complex plane of some complex parameter. Since, as seen below, the 
coalescence of just three levels has particular attractive features - for the coordinate systems 
used conventionally a distinction between left and right seems possible - the challenge to 
implement an experimental arrangement may just fall within reach of realisation. 

We recall that the wave function at a usual EP2 has - for complex symmetric matrices - a 
: ixed phase relationship of its components [14] that can be interpreted as a form of chirality 



Considering two coupled dissipative oscillators a particular EP2 specifies uniquely 



which of the two oscillators is leading by the phase tt/2. In this particular mode, that is at 
the EP2, the two oscillators thus specify an orientation in one-dimensional space by simply 
placing them on a line and using the convention that an arrow points from the oscillator 
with the leading to the one with the lagging phase is| . 

When three levels are coalescing the structure becomes much richer in comparison with 
an EP2, yet much of the structure turns out to be generalisations that could have been 
expected in hindsight. Of course, the simplest form of an operator giving rise to three levels 
coalescing is a three-dimensional matrix. Let us assume that any triple of the parameters in 
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- with U a general three dimensional orthogonal matrix parametrised by three angles -, is 
so chosen that the full problem 

Ho + A#! (2) 

has an EP3. If the parameters are all real (except for A) such EP3 will occur at complex 
conjugate values of A. Denoting such point by A c , the three levels are connected by a third 
root branch point (see appendix) and there exists the expansion: 



E j (X) = E c + ^2a k (^/X-X c 



(3) 



fc=i 



where the label j = 1,2,3 is specified by the first, second or third Riemann sheet of the 
third root in the A-plane. As a consequence, for small values of |A — A c | the three complex 
energies Ei(X) form an equilateral triangle in the energy plane. The orientation of the 
triangle depends on arg(A — A c ) and on the parameters of the specific problem ([1]), which 
determine the complex value ai in ([3]). Generically we can order the energies according 
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FIG. 1: The two basic positions of the three levels for small | A — A c | in the lower energy plane. The 
drawing is schematic in that the imaginary parts of two energies are equal: tilting either drawing by 
less than 30° still represents the specific case as discussed in the text. The lagging (120) and leading 
(-120) phases of the respective wave functions are indicated relative to the most left point (smallest 
$t(E), denoted E\ in the main text). The EP3 lies in the centre of the equilateral triangles. 

to their real parts, that is ^t(Ei) < ^.(E 2 ) < $l(E 3 ) (we dismiss the possibility that two 
real parts are equal as non-generic). There are two basic groups of orientation, depicted 
schematically in Fig.l: in (a) ^(-E^) is smaller than the imaginary parts of the other two 
energies - in other words it has the largest width - (recall that the energies all have negative 
imaginary parts), whereas in (b) Q(E 2 ) has the largest imaginary part (smallest width). 

Now we turn to the eigenf unctions. It is three eigenfunctions that become aligned when 
approaching the EP3. There are the expansions 



|^(A)> = |Vfcps) + ^(^A^X)^). 



(4) 



k=i 



Note that, as in ([3]), the label j = 1, 2, 3 is again specified by the Riemann sheet of the third 
root. For any A 7^ A c the eigenfunctions form the usual bi-orthogonal complete system, viz. 



E 



(^(A)|Vi(A)> = Nj(X)S i; j 
l^(A))(^-(A) _ _ 



r (^(a)i^-(a)) 

It can be shown (see appendix) that the scalar product §5§ vanishes as 



and similarly 



AT.(A) ~C-(A-A C )3 for A^A C 



(i>j(X)\ipEP3) ~V (A- A c )a for A — > A c 



(5) 
(6) 

(7) 
(8) 



with some constants (,rj being independent of j. 

It should be noted that the structure of the eigenvectors at an EP3 is slightly more 
involved as there are three vectors that coalesce into \1pEP3) when A — > A c . In view of the 
result ([HD the expansion (J3J) implies that not only is 

{^EP?\lpEP?>) = 



but also 



(^EP3|0l) = 



where |0i) occurs in the first order term in (jlj). 

As a consequence, and in contrast to an EP2, the eigenfunction itself does not a priory 
bear a specific chiral phase structure at an EP3. It is rather in its immediate neighbourhood 
where the chiral phase structure is revealed. Similar to the reasoning in 
of \iI>ep) in terms of the normalised basis 



llj an expansion 



lV>i(A)> 



IXi(A)> = - 

\/<%(A)|^(A)> 

yields the proper phase relation. Indeed, while it always possible for A 7^ A c to write 
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\1pEp3) = J2cj(X)xj(X) 



(9) 



(10) 



identically in A, it is in particular for A ~ A c 



2iwj 2 /3 



V 



,2iirj 3 /3 



(11) 



c 2 (A) 
\c3Wj 

with ji, j 2 , j 3 being a fixed permutation of 0, 1, 2 specified below and £ a complex constant. 
Recall that the normalised state vectors XjW blow up as l/^A — A c thus yielding a finite 
expression for \1pEP3) i n (|10|) . The result is due to the analytic structure of the third root 
branch point and is thus independent of a particular basis. As for an EP2, the fixed phase 
relations between the eigenstates in the immediate vicinity of the EP3 are of significance. 
This result is of course related to (jlj), the appealing feature of ( fTTT) lies in its similarity with 
the situation for an EP2. Note, however, that for an EP2 the phase difference is tt/2 while 
here it is 2ir/3. This means that, in contrast to the EP2, where only a fourfold loop around 
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FIG. 2: Two perspective views of the three dimensional helix when the width of E% is larger than 
the other two. The coordinates of the points denoted by are (cos($fe), sin($fc), ^t(Ek)) with 
$1,2,3 = (0°, 120°, 240°). 



the singularity restores the eigenf unctions, for the EP3 the eigenvectors are retrieved after 
three loops in the A-plane just like for the corresponding energies. 

The interesting question is now the precise association of the phases of the eigenf unctions 
with the complex values of the energies, that is with the frequencies and their widths. This 
association is illustrated in Fig.l. Choosing as reference point the eigenvector of the energy 
with smallest real part (frequency), then the eigenstate of the next frequency has a lagging 
phase of 120° if (and only if) the width is larger than the other two; the phase of the 
eigenstate with largest frequency would then be leading by 120°. In turn, if the width of 
the middle frequency is smallest, the role of leading and lagging phase is swapped among 
the two states with the larger frequencies. This result is demonstrated below in a specific 
setting and confirmed numerically in numerous general examples. 

The strict phase relations associated with the positions of the frequencies and widths 
of the three levels allows the interpretation of clear distinction between a left hand and 
a right hand helix. Fig. 2 illustrates how a particular right handed helix is generated: in 
a three dimensional coordinate system the points (cos(<&fc), sin(<E> fe ), $t(E k )) with $1,2,3 = 
(0°, 120°, 240°) invoke an oriented helix. The example refers to the case where the width of 
the middle frequency is largest, it generates a right handed helix. In turn, if the width of 
the middle frequency is smallest, the helix will be left handed. 
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FIG. 3: Suggested arrangement, arbitrarily chosen, of a microwave cavity to invoke a threefold 
coalescence of levels. As for an EP2 the opening of the gaps and some Teflon piece in one or two 
of the chambers can be used as parameters to steer the system into an EP3. To avoid systematic 
degeneracies of the levels there should be no symmetry in the geometry of the cavity. 

SPECIAL SETTING 



The all important question is: can a threefold coalescence be arranged in the laboratory 
and are the energies and phases amenable to measurement? One suggestion could be a 
setting similar in spirit to the microwave experiment for an EP2. I could imagine a setup 
as illustrated in Fig. 3. Again, this can be simulated by the simple matrix 



( ei si s3^ 



H 



si e 2 s2 
s3 s2 e 3 



(12) 



where the Sj give the couplings (gaps in Fig. 3) and the tj are the (complex) energies. To 
facilitate matters, we consider in particular the problem 



^ s3^ 



H + si ^ 



( n 1 n\ 



+ A 



e 2 si 
^s3 s2 e 3 J 

and choose the couplings s 2 and s 3 such that H has a threefold coalescence. The choice 



1 

1 




(13) 



yj u uy 



s2 



s3 



±4 



(ei - 2e 2 + e 3 ) 3 
27(ei - e 2 ) 

(-2ei + e 2 + e 3 ) 3 
27(ei - e 2 ) 



(14) 
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achieves this goal with 



E^ = ^e 1 + e 2 + e 3 ) 



The perturbation by A splits the coalescence into three levels and we obtain to lowest order 
(with A c = 0) 



/ / 2 1 /y-(-2e 1 + e 2 + - 2e 2 + e 3 ) 3 /n; n ,,x 

Ei = E r + - — — V|A - A c | exp(2«(j - l)7T/3) 



3(ex - e 2 )V3 



+ 0((A - A c ) 4 / 3 ) 



(15) 



and, using the notation (j3J), the corresponding (unnormalised) eigenvectors are 



\EP) 



\ 



\4) 



V^(ei^£2) 
^/3(ei-e 2 ) 

V 1 / 

/ -Ql/3 ^61-262 + 63 \ 

%/3( ei -e2) 5 / 6 
ol/3 V-2ei+e 2 +e 3 
v^Cei-ea) 5 / 6 



(16) 



V 
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exp(2i(j - 1)tt/3), j = 1,2,3 



(17) 



/ 



and similar algebraic expressions for the higher orders. 

The essential point here is the fact that - to lowest order - the three complex eigenvectors 
\4>{) differ only by a well defined phase; the precise specification of the particular eigenvector 
including its all important phase is given by its association with the levels Ej as indicated in 
Fig.l: the width of the second energy of the ordered frequencies characterises the orientation 
of the phases exp(2ij7r/3); if QE 2 is smallest or largest a left or right hand helix is obtained, 
respectively 



CONCLUSION 



A measurement of all three frequencies and their associated widths will thus predict the 
phase of the wave function in the immediate vicinity of the EP3. We recall that the single 
eigenstate \i^epz) does not bear such information. The relevant phases sit in the part of the 
eigenstate that is switched on when moving away from the EP3 in the A-plane. An experi- 
mental verification of this subtle behaviour would thus lead to a clear chiral characterisation 
in three dimensional space. We stress that there is no a priory handedness in a setting as 
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suggested for instance by Fig. 3, neither is there any involvement of weak interaction. There 
is, however, the direction of time that specifies the various widths of a dissipative system. 
In other words, we here suggest that the arrow of time can invoke chirality. 
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APPENDIX 

For an iV-fold coalescence of eiegnvalues the N levels are connected by a branch point 
of iV-th order. This follows from the fact that - for an iV-dimensional matrix of the form 
H\ = H + \Hi - the determinant of \H\ — IE\ vanishes linearly in the variable A while the 
requirement of iV coalescing levels entails an iV-fold vanishing in the variable E, that is the 
following set of equations is to be satisfied simultaneously 

A k 

' \H x -IE\=0, k = 0,...,N- 1. 



dE k 

This is possible only if 

oo 



m=l 



Here the N eigenvalues are given by the values upon the N sheets of the iV-th root. 
The N eigenstates coalesce likewise into one eigenvector. There is the expansion 



|V(A)) = \4>epn) + Yl l<U( 



c ) ■ 



m=l 

Note, however, that for the scalar product the following behaviour prevails 



JV-l 
JV 



mX)\^EPN)^(X-X c ) 

as follows from considering 

(E(\))-E c ){<l>(\)\if> E p N ) = (^(X)\H x -H Xc \^ EPN ) 

= (X-X^WIH^epx). 
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The right hand side vanishes linearly when A — > A c while (E(X)) — E c ) ~ y/X — A c . Note 
that these analytic properties imply the relations 

{^EPN\(pm) = for m = l,...,N-l 

and 

(0 m |0m') = for m + m'<A^-2. 
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